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Abstract 

^ Spin Foam Models are supposed to be discretised path integrals for quantum gravity constructed from the 

Plebanski-Holst action. The reason for there being several models currently under consideration is that no consensus 
has been reached for how to implement the simplicity constraints. 
' ' Indeed, none of these models strictly follows from the original path integral with commuting B fields, rather, by 

^ some non standard manipulations one always ends up with non commuting B fields and the simplicity constraints 

^ become in fact anomalous which is the source for there being several inequivalent strategies to circumvent the 

^ associated problems. 

In this article, we construct a new Euclidian Spin Foam Model which is constructed by standard methods from 
the Plebanski-Holst path integral with commuting B fields discretised on a 4D simplicial complex. The resulting 
c I ^ model differs from the current ones in several aspects, one of them being that the closure constraint needs special 

care. Only when dropping the closure constraint by hand and only in the large spin limit can the vertex amplitudes 
1— H of this model be related to those of the FK-^ Model but even then the face and edge amplitude differ. 

^ Curiously, an ad hoc non-commutative deformation of the B variables leads from our new model to the 

Barrett-Crane Model in the case of 7 = oo. 



* Unite mixte de recherche (UMR 6207) du CNRS ct des Universitcs dc Provence (Aix-Marseille I), de la Meditarranee (Aix-Marseille 
II) et du Sud (Toulon- Var); laboratoire affilie a la FRUMAM (FR 2291). 
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1 Introduction 



Loop Quantum Gravity (LQG) is an attempt to make a background independent, non-perturbative quantization of 4- 
dimensional General Relativity (GR) - for reviews, see [HElll]. It is inspired by the formulation of GR as a dynamical 
theory of connections [4|. Starting from this formulation, the kinematics of LQG is well-studied and results in a 
successful kinematical framework (see the corresponding chapters in the books [Ij), which is also unique in a certain 
sense [5]. However, the framework of the dynamics in LQG is still largely open so far. There are two main approaches 
to the dynamics of LQG, they are (1) the Operator formalism of LQG, which follows the spirit of Dirac quantization of 
constrained dynamical system, and performs a canonical quantization of GR [6l[7]; (2) the Path integral formulation 
of LQG, which is currently understood in terms of the Spin- foam Models (SFMs) [31 111 [TOl E] . The relation between 
these two approaches is well- understood in the case of 3-dimensional gravity |12j , while for 4-dimensional gravity, the 
situation is much more complicated and there are some attempts jl3) for relating these two approaches. 

The present article is concerned with the following issue in the framework of spin- foam models. The current spin- 
foam models are mostly inspired by the 4-dimensional Plebanski formulation of GR [2] (Plebanski-Holst formulation 
by including the Barbero-Immirzi parameter 7), whose action reads 

Spn[A, B,^]:= j {b + I^^bY AFij + \j A\ B'J^B^^s'^ cukl (1.1) 

where i? is a so(4)-valued 2-form field, F := dA + A/\ A\s the curvature of the so(4)-connection field A and ip°'f^'^^ = 
^p[aPVnS\ jg densitized tensor, symmetrized under interchanging [a/3] and YlS\, and traceless Cap-yS'-p'^^''^ = 0. For the 
illustrative purposes of this article, we consider only Euclidean GR in the present article, however, the lessons learnt 
will extend also to the Lorentzian theory. One can show that the equations of motion implied by the Plebanski-Holst 
action are equivalent to the Einstein equations of GR. Moreover, if we consider formally the following path integral 
partition function of the Plebanski-Holst action and perform the integral of tp"^'^^ 

Z:= J [DA DB Dif] e'^PHl^.s.^] ^ J ^5] g (eijjKLBijB^^s'' " Ve^p^s/^^-) f{B+^*B)" af,, (^3) 

we obtain the partition function of BE theory |15j whose paths are, however, constrained by 20 Simplicity Constraint 
equations 

eijjKLBijB^^g^ - -^Veapjs (1-3) 

The point of this formulation is of course that the path integral quantization of BE theory has been well-understood 
in the context of spin- foam models (see the corresponding chapters in [1 ) and thus the idea is to rely on those results 
and to implement the simplicity constraints properly into the partition function of BE theory. We remark that even 



for Euclidian gravity, the partition function (1.2) is unlikely to be derived from the canonical formulation because of 



the presence of second class constraints which affect the choice of the measure in ( 1.2 ), see the first and third reference 
in [13] for a detailed discussion. Since in current spin foam models the proper choice of measure is also regarded as 
a technical problem of higher order and as we want to draw attention to a different issue for the current spin foam 
models, we also will not deal with the measure issue in this article and leave this for future research. 

The partition function of BE theory, after discretization on a 4-dimensional simplicial complex /C and its dual 
complex K,* , can be expressed as a sum over certain spin-foam amplitudes. Here a spin-foam amplitude is obtained 
by (1) assigning an S0(4) unitary irreducible representation to each triangle / of /C (we label the representation by 
a pair (j'^ for each triangle); (2) assigning a 4-valent S0(4) intertwiner to each tetrahedron i of /C (we label the 
intertwiner by a pair {if,i^) for each tetrahedron). Then the partition function of BE theory can be written as 

ZMIC)^J2 E ndini(j;)dim(j7)n{l5j}so(4)(j/'*?) (1-4) 
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where the 15j-symbol is the 4-simplex/vertex ampUtude corresponding to the 4-simplex ct. The partition function 



Zqp turns out to be independent of the triangulation JC. Clearly, as shown explicitly in 1.2 1, in order to obtain the 
partition function for quantum gravity as a sum of spin-foam amplitudes, one has to impose the simplicity constraint 
in the BF theory measure. When doing that, the resulting partition function is no longer triangulation independent 
(as it should not be because GR is not a TQFT) and thus one should in fact consider all possible triangulations and 
not only simplicial ones. This is also necessary in order to make contact with the canonical LQG Hilbert space which 
contains all possible graphs and not only 4-valent ones. This has been recently emphasised in [8] and the current spin 
foam models already have been generalised in that respect. We believe our model also to be generalisable but will not 
deal with this aspect in the present work as this would draw attention away from our main point. 

Essentially, the very method of imposing the simplicity constraint defines the corresponding candidate spin-foam 
model for quantum gravity which why its proper implementation deserves so much attention. Currently the three most 
studied spin-foam models for quantum gravity (in Plebanski or Plebanski-Holst formulation) are the Barrett-Crane 
Model 0, the EPRL Model (TU], and FK.^ Model [TT]. These three, a priori, different models are defined by three 
different ways to impose simplicity constraint on the measure of the BF partition function Zbf- We will review these 
different methods of imposing the simplicity constraint briefly in what follows. 

First of all, in the context of the discretized path integral, the simplicity constraint also takes a discretized expres- 
sion. For each triangle / we define an so(4) Lie algebra element Bf which corresponds to the integral of the two form 
B over the triangle /. Then in terms of the Bf for each 4-simplex a the discretised simplicity constraints read 

^ijklBY B^,'" — 0, /, /' belong to the same tetrahedron t (1-5) 
eijKLB"Bf^^ ^ eijKLB"Bf^^, f„ /■ belong to the two different tetrahedrons in a (1.6) 

The Barrett-Crane Model, the EPRL Model, and the FK^ Model all explicitly impose the first type of simplicity 



constraint Eq.(1.5l, called tetrahedron constraint, in some way to the spin- foam partition function of BF theory. On 



the other hand, all of them replace the second type of simplicity constraint, called 4-simplex constraint Eq.(1.6) by 
the so called Closure Constraint 

By = for each tetrahedron t. (1.7) 

fct 

It is not difficult to see that the closure constraints together with the tetrahedron constraints imply the 4-simplex 
constraints but not vice versa. Thus, sprictly speaking, imposing the closure constraint constrains the BF measure 
more than the classical theory would precribe. It is unknown and also beyond the scope of the present paper whether 
this replacement is harmless or is in conflict with the classical theory. In this paper, as we are merely interested 
in comparing the standard way of imposing the simplicity constraints (commuting B flelds) with the non standard 
methods defining the BC, EPRL and FK models (non commuting B fields), we proceed as in those other spin foam 
models and also replace the 4-simplex constraint by the closure constraint. To distinguish these two different types of 



constraints, in what follows we use the terminology "simplicity constraint" for Eq.( 1.5 ) and "closure constraint" for Eq 



(1.7). Notice that the BC Model, EPRL Model, and FK.^ Model argue that the closure constraint is "automatically" 
implemented in their spin-foam amplitude. We will come back to this argument in a moment. Because of that 
argument, in none of these models the closure constraint is further analysed. The proper implemementation of the 
simplicty and closure constraints is one of the most active research areas in the spin foam model community and there 
are many issues that yet have to be understood |16j . 

For both the Barrett-Crane Model and EPRL Model, the strategy for imposing the simplicity constraint is the 
following: In order to take advantage of the knowledge of BF Spin Foam Model, one formally takes the delta distribution 
on the B variables out of the integral over B by a standard trick known from ordinary quantum field theories: One 
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(formally) just has to replace B by S/6F because the integrand of the B integral is of the form exp{iF ■ B). Due 
to the discretization upon which F is replaced by a holonomy around a face of the dual triangulation and B by an 
integral over a triangle of the triangulation, 5/5F can be rewritten in terms of the right invariant vector fields X 
on the copy of 5*0(4) corresponding to the given holonomy with holonomy dependent coefficients. One now argues 
that these coefficients can be replaced by their chromatic evaluation (setting the holonomy equal to unity) because 
the integration over B leads to 5{F) enforcing the measure on the space of connections to be supported on flat ones. 
Clearly, this argument is not obviously water tight because 6(6'^ /SF'^) ■ d{F) may not be supported at F = 0. In 
fact it should not be if we are interested in gravity rather than BF theory. See the chapter on spin foams in the 
second reference of [T] for more details. In any case, this way of proceeding now leads to replacing the commutative 
derivations 6/5F by the non commutative right invariant vector fields X. 

An alternative argument that has been given is the following: The kinetical boundary Hilbert space of the spin 
foam path integral should be the canonical LQG Hilbert space (restricted to the 4-valent boundary graph of the given 
simplicial triangulation) and here the B field would be quantised as 6/SA where A is the underlying connection. On 
functions of holonomies this again becomes a right invariant vector field labelled by the triangles dual (in the 3D 
sense) to the corresponding boundary edges which in turn correspond to the faces of the dual triangulation dual (in 
the 4D sense) to those triangles. The physical boundary Hilbert space should therefore be the kernel of that quantised 
boundary simplicity constraints. In order to write the corresponding spin foam model, one has to define the projector 
on that physical Hilbert space. To do this properly, one should canonically quantize Plebanski - Hoist gravity, identify 
all the first and second class constraints and define the projector via Dirac bracket and group averaging which then 
leads to a spin foam path integral. How complicated this becomes if one really performs all the necessary steps is 
outlined in [T3]. However, this is not what is done in [TD]. The first observation is that since the spin foam path integral 
naturally involves SO (4), the kinematical boundary Hilbert space is naturally also in terms of SO (4) spin network 
functions. One now studies the restrictions that the simplicity constraints impose on the spins and intertwiners of the 
boundary S0(4) Hilbert space spin network functions. The detailed structure of these restrictions suggests a natural 
one to one map with spin network states in the canonical SU(2) Hilbert space. Finally, using locality arguments, one 
conjectures that these restrictions should not only hold on the boundary but also in the bulk of the BF S0(4) spin 
foam model. See |30j for a particularly simple and clear exposition of this procedure. It has recently been criticised in 
[l6] on the ground that the BF symplectic structure and the LQG symplectic structure have wrongly been identified 
in the afore mentioned identification map. 

In any case, whether or not the map is the correct correpondence, the simplicity constraints were again quantised 
as non commuting (anomalous) constraints. If one understands the kernel in the strong operator topology then 
one obtains the BC model, if one understands it in the weak operator topology (Gupta - Bleuler procedure) one 
obtains the EPRL model. Because of the anomaly, imposing the constraint operators strongly apparently makes the 
Barrett-Crane Model lose some important information about non-degenerate quantum geometry |17| . Imposing the 
constraints weakly is less restrictive and thus may lead to a better behaved model. More in detail, first of all the 



quadratic expression of the simplicity constraint Eq.(1.5l is replaced by a linearized expression. It is given by asking 



that for each tetrahedron t, there exists a unit vector u^, such that 



*Byutj = (1.8) 



The equivalence of the linearized simplicity constraint Eq.(1.8) with original simplicity constraint Eq.(1.5l will be 
reviewed in Section [s] (in the gravitational sector of the solution). In the original construction of EPRL spin-foam 
model in [10], the unit vector is gauge fixed to be S'^'^ , and a "Master constraint" Mf :— J^j^f^f defined 
(to replace the cross-diagonal part of the simplicity constraint Eq.(1.5)), where Ci := *B°J from Eq.(1.8|. The 



corresponding "Master constraint operator" is defined by replacing Bf by right invariant derivatives. This Master 
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constraint solves the problem of non-commutativity/anomaly of the quantum simplicity to a certain extent, because a 
single Master constraint replaces all the cross-diagonal components of Eq.( 1.5 ). Moreover the diagonal part of Eq.( 1.5 1 
and this Master constraint operator restrict the Hilbert space spanned by the 4-valent SO (4) spin- networks to its 
subspace, which can be identified with 4-valent SU(2) spin-networks and thus can be imbedded into the kinematical 
Hilbert space of LQG. For each of these SU(2) spin-networks, the SU(2) unitary irreducible representations labelled 
hy k G has the following relation with the original S0(4) representations on all the boundary edges dual to the 
boundary triangles 

Here the Barbero-Immirzi parameter 7 can only take discrete values, i.e. 

If|7|>l: .-'^-"'^ 



If|7|<l: 7=^^-4 (1.10) 
More importantly, the recent results in [18 show that the boundary Hilbert space used in the EPRL Model solves the 



linear version of simplicity constraint Eq.( |1.8[ ) (and the closure constraint Eq.(1.7)) weakly, i.e. the matrix elements 
(with respect to the boundary S0(4) Hilbert space) of the constraint operators vanish on the space of solutions 

/, C*/'^ — 0, for all /, /' in the Hilbert space of solutions. (l-H) 

in contrast to the strong implementation of the constraints in the Barrett-Crane Model. Finally the (Euclidean) EPRL 
spin-foam partition function is expressed by 

w(/c) = E E ndi-(^/)nE{i5j}so(4) '*^) n (^t (i-i2) 

{fc/}/{it}t i " if 

where for each spin- foam amplitude, an SU(2) unitary irreducible representation fc/ is assigned to each triangle /, 
satisfying the relation Eq.(1.9), and an SU(2) 4-valent intertwiner it is assigned to each tetrahedron t. Here 

E{i5j}so(4) (j/ n q^ (^/ (1-13) 

is the 4-simplex/vertex amplitude for the EPRL Model, where /*| are a fusion coefficients defined in [TU] . 

The FK^ Model follows a different strategy to impose the simplicity constraint, namely by using the coherent 
states for SU(2) group [111 [20]. Given a unitary irreducible representation space of SU(2), the coherent state is 
defined by 

i 

:=n|.?,j) = E |j,m)7r^^^.(n) n G SU(2) (1.14) 

We then immediately have the resolution of identity on 

lj=dim(j)/ dn|j,n)(j,n| (1.15) 

J SU{2) 

This coherent state has a certain geometrical interpretation, which can be seen by computing the expectation value 
of the su(2) generator {pi are Pauli matrices) 

(j, n\ X \i, n) = (j, n\ > | j, n) Oi = jna^n^^ (1.16) 
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If we identify the Lie algebra su(2) with M.^, we can see that the coherent state \j,n) describes a vector in M.^ with 
length j, its direction is determined by the action of n on a unit reference vector (the direction of a^). From the 
expression na^n^^ we see that n can be parameterized by the coset SU(2)/U(1) — S^. In addition, the integral in the 
resolution of identity is essentially over SU(2)/U(1) = S^. It is not hard to show that the (Euclidean) BF partition 
function can be expressed in terms of the coherent states (we write {g'^,g~) for each S0(4) element, for an 

SO (4) unitary irreducible representation) 



(tj) (a/) 



9ta9tt 



9t(y9crt 



Jj,n;,,) (1.17) 



where [gtt^9sigt) ^ S0(4) holonomy along the edge from the center of 4-simplex a to the center of tetrahedron t. 



Then the strategy of imposing simplicity constraint in FK-y Model is to use the interpretation (1.16) of the coherent 



state labels j^n^r3(n^) ^ as the self-dual/anti-self-dual part X^^ of the so(4) variable Btf associated with a triangle 
/ seen from a tetrahedron t. (More precisely, we know that the previously defined Bf can be decomposed into self-dual 



The interpretations of j^nfjT3{rifj) ^, Xfj are considered as the parallel transport of 



and anti-self-dual part Xj^ 

X^ from the center of triangle / to the center of tetrahedron t, i.e. X^j: = gf^j^X^g^^, where is the holonomy along 
the edge from the center of triangle / to the center of tetrahedron t). That is, the simplicity constraint is imposed on 
the coherent state labels, which results in the following restrictions: 



7 + 1 



7 



1 



and 



i^tf^^tf) = {ntfh^,j,utntfh^l^), 
{r4f,n-j) = {ntfh^,f,utntfh^^^e), 



for -1 < 7 < 1; 

for 7<— lor7>l. 



(1.18) 



where Ut is some normal to t, h^^j: takes values in the U(l) subgroup of SU(2) generated by CT3 and e — ia2- In more 
detaiul, the proposal is then to simply replace in (1.17) by these expressions and the Haar measure dn'^j.- dri^j 



by the Haar measure dritf dut dh^^^. We emphasize that this is an interesting but non standard procedure: While 



the identification of the coherent state labels 



with the so(4) variables Btf is certainly well motivated, the 



resulting expression does not arise by integrating out the B fields in the presence of the delta distributions enforcing 



the simplicit constraints. Rather, in (1.171 the B fields have already been integrated out. To restrict measure and 



integrand by hand afterwards according to (1.18) is not obviously equivalent with the standard procedure of solving 



the 5— distributions. One would hope that the resulting procedures coincide in the semiclassical or the "large-j" limit 
[21] . Indeed, the "large-j" limit result in Section |4] will support this expectation. Finally the spin-foam partition 
function of FK.^ Model coincides (at least up to a slight change of edge amplitude) with EPRL partition function 
when the Barbero-Immirzi parameter —1 < 7 < 1. However when 7<— lor7>l, FK-y partition function is rather 
different from the EPRL partition function. Here we only show explicitly the 4-simplex/ vertex amplitude of FK-y 
model when 7<— lor7>l 



(1.19) 



(cr,t) 



Here although the relation between 



7 + 1 



7-1 



(1.20) 



is the same as in EPRL Model, in FK..y model for 7 < —1 or 7 > 1, there are some additional degrees of freedom 
associated with the label ktf, which are the values of spins from the coupling of J/ j ^t/ could take values 
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in \j^ — jf \r ' ' +jf ■ The final partition function is obtained by summing over jj , it, and ktf with, some measure 
factors (see [Tl] for details). 



In the previous three paragraphs, we briefly revisited the main strategies of imposing simplicity constraint in 
Barrett-Crane, EPRL and FK^ Models. We have seen that these in general different spin-foam models came from 
two different ways of imposing simplicity constraint, i.e. Barrett-Crane and EPRL Model quantize the simplicity 
constraint as operators and imposed them (strongly or weakly) on the boundary spin-networks, while FK^ Model 
imposes the constraint on the coherent state labels. However, as we have reviewd, none of the three models is 



derived from the original path integral formula Eq.( 1.2 1 of the Plebanski action (or the discretized version of the path 



integral) without using some non standard methods. Therefore a natural question arises: 



Is any of those three spin-foam models consistent with the path integral formula Eq.(1.2) and its discretized version? 
This question is non trivial because in all three types of models one deals with non commutative B fields and 
simplicity constraints as operators on some Hilbert space while the original path integral is in terms of commutative 
c-number variables so that anomalies cannot arise. Because of this issue, it is interesting to investigate what kind 
of spin-foam model we will obtain, if we start from the (discretization of) the path integral formula Eq.( |1.2| with 
commutative B^' variables. It is also interesting to find some possible bridges linking the (discretization of) the path 
integral formula Eq.(1.2) with commutative B^"^ variables to the existing spin-foam models using non-commutative 
B^"^ variables. 



In this article, we consider the discretization of the path integral formula Eq.(1.2|, which will be Eq.(2.1). As 
announced in [21], in contrast to the Barrett-Crane, EPRL, and FK.^ Models, we always consider the variables B^^ 
as commutative c-numbers. The simplicity constraint (and closure constraint) is (are) imposed by the c-number delta 
functions inserted in the path integral formula, which one gets by integrating over the Lagrange multiplier and which 
constrain the path integral measure. In our concrete analysis in Section |4j the most important difference between our 
derivation and the derivation in any of Barrett- Crane, EPRL, and FK-^ Models is the following: in any of Barrett- 
Crane, EPRL, and FK-y Models, one always imposes the respective version of the simplicity constraint constraint 



on the BF spin-foam partition function Eq.(1.4) or (1.17) after integration over B^"'. This feature is essentially the 



reason why it is difficult to find a relation between the simplicity constraint imposed in any of Barrett-Crane, EPRL, 



and FK^ Models and the simplicity constraint in the path integral formula Eq.(1.2). By contrast, our derivation in 
Section [4] will not start from the spin- foam partition function of BF theory, but instead we impose the delta function 
of the simplicity constraint (and closure constraint) before the integration over B", and we will see that solving these 
constraints gives rise to a non trivial modification of the path integral measure. 

As also announced in regarding the B^"^ variables as commutative c-numbers also makes the treatment of 



closure constraint different. We know that the closure constraint Eq.(1.7l is necessary in order that the full set of 



simplicity constraint Eq.(1.5) and (1.6) is satisfied. In Barrett-Crane Model the closure constraint is argued to be 
automatically satisfied by the S0(4) gauge invariance of the vertex amplitude. However, as shown in [31], this is only 
true after performing the Haar measure integrals which essentially project everything on the gauge invariant sector. 
It is clear that the closure constraint must be imposed before performing the integral over the connections. In the 
EPRL Model, the argument is improved in that both simplicity constraint and closure constraint vanishes weakly on 
the EPRL boundary Hilbert space pLBj. Moreover, in [22], it is shown that in both EPRL and FK-y Model, the closure 
constraint can be implemented in terms of geometric quantization and by the commutativity of the quantization and 
phase space reduction [53]. As defined, an additional closure constraint would be redundant for both EPRL and FK.^ 
Model, since they are already on the constraint surface of closure constraint (if one interpret the coherent state labels 
to be the B^'^ variables), although the original definitions of both models didn't impose closure constraint explicitly. 
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We feel that this is again due to the fact that the Haar integrals have already been performed. In our analysis we find 
that the implementation of closure constraint gives non-trivial restrictions on the measure. 

In order to understand what happens when one ignores the clsoure constraint and to follow more closely the 
procedure followed by existing spin foam models, in section [4j we first consider a simplified partition function 
■Z^Simpiified(A^) in which the delta functions of closure constraint is dropped, and derive an expression of 2'sii„piified(A^) 
as a sum of all possible spin-foam amplitudes (constrained only by the simplicity constraints) . Then we compare with 
the full partition function Z{1C) with the delta functions for closure constraint implemented, and we find the delta 
functions of closure constraint non-trivially affect the spin-foam expression of partition function. Surprisingly, all the 
spin- foams (transition channels) admitted in the simplified partition function .^simplified (^) also contribute to the full 
partition function Z{IC) (with some changes for the triangle/face amplitude and tetrahedron/edge amplitude). 

Another key feature of our derivation is a different discretization of the BF action. Here we first break the 
faces dual to the triangles into wedges (see FIG{T]) and then write the discretized BF action in terms of the holonomies 
along the boundary of the wedges. Here, as usual, a wedge in the dual face / is determined by a dual vertex or 
original 4-simplex a and thus denoted by (cr, /). Its boundary consists of four segments defined as follows: The 
original (piecewise linear) 4-simplex has a barycentre a which is the dual vertex. The dual edges connect these 
barycentres. A pair of dual edges e, e' adjacent to the same dual vertex defines a face. Conversely, given a face and a 
dual vertex which is one of the corners of the face, we obtain two dual edges. These are dual to two tetrahedra t, t' of 
the original complex. The boundary of the wedge {a, f) is now given by (it, e) o (e, /) o (/, e') o (e', a) where the hat 
denotes the respective barycentres. In an unfortunate abuse of notation which exploits the duality one also writes 
this as (cr, t) o (t, /) o (/, t') o (t'j a). Using this notation we have (cf. FIGjl} 




(1.21) 



where Fi^a.f) is the curvature 2-form integrated on the wedge determined by (cr, /) and i,t' respectively are the afore 
mentioned unique tetrahedra (or dual edges). This starting point results in some new structures in the resulting 
spin-foams: 

• In contrast to the existing spin-foam models, where the S0(4) representations {jj^ were labeling the faces /, 
the new spin-foam model derived in Section |4] have S0(4) representations (j'^pj'^j:) labeling the wedges, i.e. a 
dual face / having n vertices (corners) in general has n different pairs (j'^pj'^j), one for each wedge determined 
by the vertex dual to a. However in the large-j limit, the triangle/ face amplitude is concentrated on S0(4) 
representations — f) for any vertices cr, cr' of the same face /. 

• Two neighboring wedges (cr, /) and (cr', /) of a face / share a segment {t, /) (cf. FIGjl]) whose end points are the 
center of the face / and the center of the edge dual to the tetrahedron t = cr n cr'. For each segment (i, /) there 
is an SU(2) representation ktf "mediating" the S0(4) representations on the two neighboring wedges, {j^^j'^) 
and {j'^,f,j'^,f), in the sense that ktf has to lie in the range of the joint Clebsh - Gordan decomposition of 
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jtf^jaf ^ndj+f®j^,f (ci. FIGQ, thus 

ktfe[\j+f-j;j\r-- ,j+/+j7/}n{|j;,;-j7'/l'--- ^J^'z+V/}- (i-22) 

It is interesting to see the relations between the new spin-foam model derived here and the existing spin-foam 
models e.g. Barrett-Crane, EPRL, and FK.^ Models. From the analysis in Section |4j we find that, firstly, in the 
large-j limit the spin-foams in our new model .^simplified (A^) reduces to the spin-foams in FK^ Model (with identical 4- 
simplex/ vertex amplitude but different tetrahedron/edge and triangle/face amplitudes) at least for I7I > 1. Secondly, 
in Appendix [Aj we study the non-commutative deformation of the partition function Eq.(2.1), in order to study 



how the non-commutative nature of the B^'^ variables in the existing spin-foam models emerges in our commutative 
context. The non-commutative deformation we employ here comes from a generalized Fourier transformation on the 
compact group [21] (the deformed partition function will be denote by Z^(/C)). With this deformation, we find that 
the closure constraint really becomes redundant when we set the deformation parameter a — £p, while the redundancy 
is hard to be shown with a general deformation parameter. With the setting of the deformation parameter a = i^, 
we show that the non-commutative deformation of our new spin-foam model leads to Barrett-Crane model when the 
Barbero-Immirzi parameter 7 = 00. This result explains how the non- commutative nature of the B^"^ variables in 
Barrett-Crane model relates to the commutative context of our new spin- foam model in Section |4| and also explains 
to some extent the reason why in the Barrett-Crane model the closure constraint is redundant (such an explanation 
also appears in the first reference of [IS] from the group field theory perspective). On the other hand, the relation with 
EPRL Model and FK.^ (I7I < 1) is still veiled. What we know is that the allowed spin-foams (transition channels) in 
EPRL Model form a subset of those allowed in our new spin- foam model (with the same 4-simplex/ vertex amplitude 
but different different tetrahedron/edge and triangle/face amplitudes) and this fact also holds for FK.^ Model for any 
7. All above relations between various spin-foam models are summarized in the following diagram, where the sets 
{Z...} are the collections of spin- foams (transition channels) which respectively contribute their partition functions 
Z...(/C): 



{ZePRl},{ZfK^,} C {^Simplified} C {Z} {Zbc} 

4, largo- J liirrit, 1^1 > 1 

{^FK|.,l>i} 

where C means the inclusion in terms of contributing spin-foam amplitudes. We will discuss the details in Section 4.2. 



2 Starting Point of the New Model 
2.1 The Partition Function 

In the last section we reviewed the approaches of simplicity constraint and closure constraint in the existing spin-foam 
models, and summarized the approach and main results of the present article. In this section, we present the detailed 
construction and analysis of our new spin-foam model. We take a simplicial complex JC of the 4-dimensional manifold 
M, where we denote the simplices by cr, the tetrahedra by t and the triangles by /. And we take the following 
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discretized partition function as the staring point for constructing the spin-foam model: 

Z(/C) := / n d^^;d^^7 n d5+.d5-, J] ^^atf^dtf 11 ^ (^*/ ' ^*/' " ^tf ' ^r/') R < E ^*/) 
/ (<y,t) (tj) t-JJ'dt t /ct 

We explain the meaning of the variables appearing in the above definition: 



(2.1) 



, Xj S su(2) are respectively the self-dual and anti-self-dual part of the so (4) flux variable By , which is the 
so(4)-valued 2-form field smeared on the triangle dual to / while 



So given two tetrahedra t,t' sharing a face /, the relation between Xtf and Xfi f is thus 

^t'f ■- 9t't^tf9u' 



(2.2) 



(2.3) 



where g^,^ = g^, ^g^^ and g^^ — {g^t) ■ This means that the "shape-matching" condition [26^ is imposed in our 
model 

Ag is the Haar measure on SU(2). g'^^^ e SU(2) is the self- dual and anti-self-dual part of the S0(4) holonomy 
along the half edge {a, t] outgoing from the vertex a while g^j:, g^^are respectively the self-dual and anti-self-dual 
part of the SO (4) holonomy along the segments (t, /) (see FIgEI). 




Figure 1: A face dual to the triangle /. The vertices of the face are dual to the 4-simplices a. Each edge of the face is 
dual to a tetrahedron t. The fundamental region (in gray) in the face determined by a 4-simplex and two tetrahedra 
is called a wedge. Each tetrahedron is shared by two 4-simplices. A triangle is shared by n simplices, where n is the 
number of vertices of the dual face. 



The delta function S (^X^j: ■ X^j, — X^^ ■ X^^^ imposes the simplicity constraint for each tetrahedron: 

^ijklB^ B^jf" — f,f' belonging to the same tetrahedron (2.4) 

while the delta function ^(^J^fct^tf^ imposes the self-dual closure constraint for each tetrahedron. Note that 
there is no closure constraint for XJ because the closure of XJ is implied by the self-dual closure constraint 
and the simplicity constraint as we will demonstrate shortly. So including it would be equivalent to multiplying 
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the partition function with a divergent constant which drops out in expectation vahies. In addition, the closure 
constraint and simpUcity constraint Eq.(1.5) imply the 4-simplex constraints fc, / € {1,2,3,4,5}): 



fijiace dual to the triangle ti D t,-, where are the 5 tetrahedra of a 



(2.5) 



Here Xfj: — g^^Xf'^gf'^ and B„f = X'^^ + X^j. In the continuum limit of Eqs.(2.4) and (2.5 1, in which the 
holonomies can be replaced by the group unit, we recover the Plebanski simplicity constraints (20 equations): 

eijKLBijB^g'^ = Ve„^^A-/4! (2.6) 

where V :— f-ij KhBlfpB^g" is the 4- dimensional volume element. Note that there are essentially 20 

constraint equations while the trace part of Eq.(2.6 1 is an identity. The solutions of the simplicity constraints is 
well- known: given a non-degenerate co-tetrad e^, there are five sectors of solutions of the simplicity constraints 

m 



i± 
ii± 

Deg 



B'' = ±e' A e"' 

B^' = ±\e"^,e^^e^ 

B'' = 



(2.7) 



The exponentials in Y[ 



^^(l+i)tr(x;3+g+g+,3+^.) 



n(., 



come from the exponential 



of the BF action, regularized in terms of wedge holonomies gatQwdva- more detail. 



M 



1 

B + -*B 

7 



^F,.^ = 



M 



1 



tr ( x+f; 



tr X+ AF^ 



tr[XJF^ 



M 



1 - 



tr {X- A P- 



tr X 



9ft9ta9at' 



tr X 



9ft9ta9at'9t'f 



where F(^a-.f) is the curvature 2-form integrated on the wedge determined by (tr, /) 
Finally we note that under the SO (4) gauge transformations: 



(2., 



± v±(u±\-l 



Xf ^ hJXJ{hj _ 



Xff ^ hfxfAhf)-^ 



(2.9) 



where /i : E — > 5*0(4); x i-> h{x) denotes a gauge transformation and ■— h{a), ht := h{t), hf{f) with a the 
baryonic centre of a etc. 

Hence the traces of the exponentials 



tr 



(^Xfgf^gt,gf^,g^j 



(2.10) 



and the simplicity constraint 



(2.11) 
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are invariant quantities while the closure constraint transforms covariantly 

Y^xtf^h, (E^t/I ^r'- (2.12) 

/ct \/ct / 

The desire to maintain gauge (co)invariance of action and constraints in the discretisation motivated to introduce 
the quantities X^^ and X^^ which in the continuum limit reduce to X^ to leading order in the discretisation 
regulator. 

• One may wonder why we do not include 5 functions enforcing the closure constraint for the "minus" sector. As 
we will see, the measure is supported on configurations satisfying X^^ — utX^jU^^ for some ut €SU(2). Thus 

E^*7 = -"t[E^t/H"' (2-13) 

/ct /ct 

is already implied by the "Plus" sector. So we could include it but that would result in an infinite constant 6{Q) 
which drops out in correlators. We assume to have done this already. 

Remark: 

It appears awkward, that here are more holonomies than B fields, suggesting a mismatch in the number of B and 
A degrees of freedom in contrast to the classical theory. Here we remark that the natural definition of the dual of a 
triangle really is the gluing of wedges (see e.g. the second reference of [T] in the notation used here and references 
therein). The boundary df is naturally a composition of the half edges [t, tr] where the hat denotes the barycentre 
of tetrahedron and 4-simplex respectively. Thus, if we would discretize the action using the holonomy around the df 
rather than around the wedges, the discretized action only would depend on the edges e — [aa D a'] D [ana', a'] and 
the properties of the Haar measure ensure that the inegrals over g^^g^^ reduce to the integrals over gf. Thus, what 
we are doing here is to approximate tr(_B/ ■ gf) by X^^e/ t^(-^/ ' 9f,c) where gf^a = 9 ft gta gat' gt' / is the corresponding 
wedge holonomy a/ier having introduced the redundant variables gta, gtf- We are aware that this presents a further 
modification of the model but it should be a mild one because both discretised actions have the same continuum limit. 
In fact we will see that in the semiclassical (large-j) limit the representations on the wedges essentially coincide so 
that effectively only the face holonomies are of relevance. It is certainly possible to define the commutative B field 
model without this step, however, it is very helpful to do so as it facilitates the solution to otherwise cumbersome 
bookkeeping problems. We leave the definition of the model without a priori introduction of wedges for future work. 

2.2 Expansion of The Exponentials 

For the preparation of the integration of the holonomies g^f amd g^, we would like to expand the factors 
g*(i±T,)tr('^/ S/t9t,T9„t/9t//) terms of the SU(2) unitary irreducible representation matrix elements 7r^„(5). So we 
define the matrix K-i^^(Y), Y e su(2), such that 



= J2 Kir^niYXnig) (2-14) 



while the expression of K^^{Y) can be obtained by 

^^fmniY) ^Jdg e^'^'-oKLnig) ^Jdg e'*'-(^^)<„(5-i) (2.15) 
Since iY = iy ■ t = y ■ a (aj are Pauli matrices, Tj = —iaj), we have the following relation 

= M,y . 0^ ^ Mg^f!'- (2.16) 
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Therefore 



dim(j) 



-x^jy) = / dg 



(2.17) 



where in the last step we made a translation g ^ ge ^S* '^. Moreover we can expand the function e'? ''''^-' by the 
SU(2) characters 



-i|!7|tr(g) 



(2.18) 



fceN/2 



Then 



fceN/2 



E /?fc(lyl)E<<(^^' 

fceN/2 i 



d.9 '^iniig)xkig) 



E /^fc(iyi)E<<( 



e'2 



fceN/2 
1 



dini(j) 



dim(jr'' 

Then plugging this result back into Eq.( |2.14 ) yields 



(2.19) 



(2.20) 



by using this identity, we have {X^ = ■ f = X^ ■ (—«(?)) 

^(l±i)tr(x±s±3±9^,,^^,) ^ 



1 

1 ± - 

7 



Xl 



Inserting this result into the expression of the partition function, we obtain 



(2.21) 



/ (^T.t) (t,/) t;/,/'ct t /ct 



E n h 
E n ^ 



/ l) ^"^j-:, (*^/ • ^9Ut^9tt'9t'f 
) (*^7 • '^9]t9ia9~t'9rf 



(2.22) 



3 Implementation of Simplicity Constraint 
3.1 Linearizing the Simplicity Constraint 

In order to implement the simplicity constraints via the delta functions 5 [x^ ■ X^, — XJ ■ XJ,^ for each tetrahedron 
it proves convenient to pass from this quadratic expression to an integral of linear expressions directly at the level 
of measures (in the gravitational sector II±). In this subsection we are dealing with a single tetrahedron t, thus we 
ignore the t label of X^j^. 



Consider the four flux variables Xf (/ = I,-- - ,4) associated with a tetrahedron t. Define the symmetric 
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matrix Z^^, := ■ X^, 1 </,/'< 4. Then Z^^., determines the X^ up to an 0(3) matrix O. Denote by L the range 
of the map {X^}-j^^-^ i-)- {-^/}/}i</</'<4 (as a subset of M.^^, L is constrained in particular by the Cauchy-Schwarz 
inequaUty). Then we can define a map Y : 0(3) (g) L K^^, {g,l) {gXi{l),gX2{l), 9X3(1), gX4{l)) where Xf{l) is 
any solution oi Iff = Xf ■ Xf. 

In the following result we drop the ± for convenience. 

Lemma 3.1. We have det((Z//') = 0. Given F : R^"^ define F : 0(3) x L ^ R by F := f oY . Then 

I f[d'XfF = [ dg [ d'°l 6 (det((Z//0)) F (3-1) 

>/Mi2 t_, Jo(3) JrxL 



/=1 



where dg is the SU(2) Haar measure (up to normalisation) and F is trivially extended off the surface det{l) = 0. 

Proof: Up to measure zero sets, Xi,X2.,X3 will be linearly independent and define a 3 metric lab = Xa ■ X^. 
Accordingly 

X4 = l'^''(Xf, ■ X4)Xa = 1°'^ lb4 Xa (3-2) 

is a linear combination of these vectors and l°'''lcb = ^h- We obtain the constraint 



'44 



■ X4 ■ X4 — l^^Uahb (3-3) 



among the Iff. On the other hand 



det((Z//0 = det ( ] = I44 det{lab - hahb/ki) = I44 [det(X^)]2 det(Z44% - Z4i^4i) (3.4) 

\ 'a4 lab J 

with l4i = X^l4a and Xf is the inverse of X^. The computation of the remaining determinant is elementary and yields 

det((Z//0) = det((Z„b))[X44 - /''''Z4a/46] (3.5) 

which is proportional to the constraint, hence det((Z///)) = 0. 

In order to write an integral over Xi, .., X4 in terms of the independent coordinates lab, Ua, ct where a parametrises 
the rotation g, we must compute the Jacobian 

T I J .fd{Xi,X2-X2,X4) 

^='^^^( dilab,a,l4a) ^' ^'-'^ 

Since only X4 depends on this immediately simplifies to 

To compute the remaining determinant we choose for instance the following parametrisation 

Xi = -y=b3 + \Jhi - Ih/hsicosi'j + x)bi + sin(7 + x)b2) 
X2 = 



ha 

I23 



bs + \Jh2 - ^23A33(cos(x)6i + sin(x)62 



V^33 

^3 = Vhsbs (3.8) 
with the Euler angles a = (0, 0,x), X € [0, 27r], 9 € [0, tt] and the orthonormal right oriented basis 

63 = (sin(^) cos(0), sin(^) sin(0), cos(^)), 61 = 63,9, 62 = 63,^/ sin(^) (3.9) 
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together with 

I X ^12^33 — ^13^23 . / N ■\/det(/ab)/33 /0 1n^ 

cos(7) = , sin(7) = (3.10) 

V (^11^33 — ^13)('22^33 — ^23) V (^11^33 ^ 'i3)(^22^33 ~ ^23) 

This defines the map Y above and the reader may check that the relations • Xi, = lab arc satisfied for any a. The 
computation of the Jacobian is much simplified by noticing that the matrix 

d{X,,X2,X,) ^^^^^^ 

9{l337 ^7 01 Xi ^22, ^23; ^11; ^12, ^13) 

consists of 3x3 blocks and is upper block trigonal with non singular matrices as diagonal block entries. Accordingly 
its determinant is the product of the determinants of the diagonal block matrices and yields after a short comutation 



the value sin(0)/(8y^det(I^). Due to the absolute value the Jacobian is thus given by 

J=^?Tr^ (3-12) 

8det(Zab) 

It is not difficult to check that for the Euler angle parametrisation we have up to a normalisation constant the following 
expression for the Haar measure 

dg = dx d(j) de sm{e)/8 (3.13) 

We can therefore finish the proof by 

d^Xi d^X2 d^Xs d^Xi F 

I dg Yl dlab n dha , 
o(3) Jb „i,^<, t\ det((U)) 

4 



/ dg f r\ dlab Y\ dhf -J^— Sihi-l'^'hakb) 

Jo{3) JbxR ^<j<3 fj-^ det((ia6j) 

f dg f n dlff, sniff,)) F 

JO(3) JBxR f<f'<4 



(3.14) 



□ 



As usual in path integrals we will not worry about normalisation constants as they drop out in correlators. The 
preceding lemma is crucial for establishing the following result. 

Lemma 3.2. For each tetrahedron t (ut G 50(^5^^ 

W 5 {x+ ■ X+ - XJ ■ XJ^ ^ S (det(X+ ■Xp'j / dut Y[ ^ (^X+ + utXJut^^ (3.15) 

/,/' = ! /=1 

in the solution sector II± of the simplicity constraint ill}/ . 

Proof: Essentially we need to prove that for all continuous function f{X^,XJ) f — 1, • • • ,4 vanishing in the 
topological solution sector I± of the simplicity constraint 

/ n d^^;d^^7 n ^ (^/ • ^7 ■ ^/") / (^/ ' ^7) 

/=1 /,/'=! 

» 4 „ 4 

= / n d'^/d^X/J (det(X+ • X+)) dutY[S (X+ + UtXjui^^ f [x+,XJ^ (3.16) 
/=i /=i 

^SO(3) is considered as the upper hemisphere of SU(2), while their Haar measure is different by a factor of 2. 
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From the left hand side, by using Lemma 3.1 we transform the coordinates from to and l^fi, constrained by 
det{lff,) = 



f=i fj'=i 

= I du+dV n d^;,,d/7,,J (det(/;,,)) -5 (det(/7,,)) R ^ (^//' " ^7/') / (^//" ^7/" 
= / du+Aui n d;;/,'5 (det(?;^,)) 5 (det(?+^,)) / "t") 

= / n d'^/ / d^r^ (det(x; • x+)) / 

» 4 „ 4 

= / n d'X+d3X7<5 (dct(X; • / d^/t [] ,5 (X+ + utXJui"^) f [X+,XJ^ (3.17) 

where we restrict ourself in the gravitational sector J/±. □ 

Notice that strictly speaking we should be using the Haar measure duf on 0(3) rather than S0(3) which is just 
the sum of two Haar measures on S0(3) twisted by a reflection so that we actually get an integral over SO{3) of a 
sum of S distributions 6{X — UtX'^u^^) + 6{X + UtX'^u^^) with u-t e 50(3). This is expected because the simplicity 
constraints do not select either of the two sectors (gravitational and topological). As usual in spin foams, we consider 
a restriction of the model to the purely gravitational sector in the above lemma. 

Here we note that the singular factor 6{det{X^j: ■ X^^',)) is essentially a 5(0) and can be divided out by an appropriate 
Faddeev-Popov procedure [11]. And the linearized simplicity constraint S (^X^ + utXJu^^^ has clear geometrical 
interpretation that for each tetrahedron t, there exists a unit 4-vector rif = {nl , , nf , n^) corresponding to the SU(2) 
element 



= / 3 ■ 4A 1 -2 (^-^^^ 



such that *B f rit / = 



Thus the constrained measure of the flux variables in Eq.(2.22) is written as (we denote gt = ut in what follows) 



\{d'X+A^Xj\{ [ Au^\{6(x++u,X;^ui^)\{6[Y.X+) (3.19) 



f t " f<Zt t fct 

can be considered as the measure f^'^A'^j 



Note that the measure (fiXf can be considered as the measure <PX^f constrained by the "shape-matching conditions" 

5(a±-5,±,a±(5±)-; 

In particular we see, that it is possible to justify the passing between the quadratic simplicity constraints employed 
by the BC model and the lineraised simplicity constraints of the EPRL and FK models respectively, directly at the 
level of measures in terms of the commuting B variables. 

3.2 Imposing the Simplicity Constraint 

In what follows we make the ad hoc restriction to the gravitational sector as mentioned at the end of the previous 
subsection. 
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Performing a polar decomposition of the variables and ^^j, we introduce the new variables G and 
Nf e SU(2) 

Xf^pfNfMNf)-' X± =p±7V±r3(iV±)-i 7V±^g±iV± (3.20) 



where pf = \\Xf\\, xf ■ f ^ Nfr^iNf)-'^ and the same for X^^. Note that given X± e su(2), iV± e SU(2) is 
determined up to a U(l) rotation /i0 S U(l), which leaves ra invariant. 



e'"^ 
e-*-^ 



(3.21) 



The associated equivalence relation is called the Hopf fibration of SU{2) = S*"^ as a C/(l) bundle over the coset 
space SU{2)/U{1) = 5^. It is convenient for given unit vector n{d,(j>) — (sin(0) cos(0), sin(0) sin(0), cos(0) to fix the 
representative N = ie^ {6,(t))crj with the unit vector e{6,(t)) — (sin(26') cos(0), sin(26') sin((/)), cos(26') parametrising a 
point on S^. 

The the linearized simplicity constraint X^j = —utX^jU^^ implies that there exists a h^^^ e U(l) for each pair of 
{t, /) such that 

p+ = p- = pf and [n+,N-^) = (Ntfh^,^ , utNtjh-^le) (3.22) 

where the diagonal U(l) invariance is absorbed into the definition of Ntf, we only take care of the anti-diagonal one 
by introducing (ptf, and 



1 

-1 



(3.23) 



We now reexpress the constrained measure in terms of the new variables p^ and . The Lebesgue measure d'^X 
can be expressed in the spherical coordinates (for any function of X) 

d^X = p^dpd^Vt p2 dp dN (3.24) 

where d^l7 is the round measure on and dA^ is the Haar measure on SU(2). 

Lemma 3.3. For any continuous function f{N^,N^) on SU(2)x SU(2), up to overall constant factor 

dA^+diV" Ss2(n-T3{N-)-^ +uN+T3{N+)-^u-^) f{N+,N-) 

dN+dN- f dN f d^ 6su(2){N+,Nh^) Ssu(2)(N-,uNh^h) f{N+,N-) {3.25) 



SU(2)xSU(2) 



2-K 

\N I 

SU(2)xSU(2) JSU(2) 

which gives 

r2iT 



Ss2[n-T3{N-)-^+uN+ts{N+)-\-^) = d<l>6su(2){N-,uN+h-^e) (3.26) 



Proof: On the right hand side of Eq.(3.25), 

r2n 



dN / deb 5su(2) [N+,Nh^') 5su(2) [N-,uNh^h) f{N+,N-) 



f(^N+,uN+h-^e) (3.27) 



dN+dN- 

SU(2)xSU(2) JSU(2) 

dN+dN- / d<j>Ssvi2){N-,uN+h24,e) f{N+,N-) 

SU(2)xSU(2) Jo 

dN+ I ( 

SU(2) JQ 
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On the left hand side, we can express the Haar measure dA^ in terms of Euler angles 



dN- ■ ■ ■ 



SU(2) 



167r2 7o 











02 / d6' sin 6* 


/ 




"'0 





(3.28) 



And the delta function u is the delta function on S"^ , which is coordinatized 

by G [0,7r] and (/)2 € [0,27r]. By explicit computation 



2ir i^TT 



de sine 5s^{N^r3,{N~y^ +uN+Ts{N+)-^u-^^ f{N^ 



Therefore the left hand side of Eq.(3.25) reduces to 



SU(2) 



47T 



dN+ / d(/) 



f[N+,uN-h^h). 



which is identical to the right hand side Eq.(3.27). 
□ 



Using this we rewrite the constrained measure up to an unimportant overall constant as 



fct 



t fct 

I.2TT 



ndp+d7v;(p+)'dp7d7V7n / d-*n^(4-^7) / '^#*/5K,"*^^t>2-,v^) 



t 'fct 

We insert this result into the partition function 



z(/c) = / n '^a^tdg-t n ^41.9*7 n M'^^t {p^) ' '^pj'^^i n n ( 

(<T,t) (tj) f t fct 

2ir 



Pf - Pf 







dq^,fS [Nr^,u,N+h~l^e) n '5( E P^K^Ky' 

t fct 

p^) tr.:, (^^;-3 (iv;)"%+ff+g+,5^ 

(*^7^3 (^7) 9jt9ta9~t'9;:f^ 



X E n 

X E n z^.;, 



Performing a translation of the Haar measure dg^ 



d5±^d(g±A^±)=d7V± 



(3.29) 



(3.30) 



(3.31) 



(3.32) 



(3.33) 



(notice that diV^^ and dA''^ are Haar measures on SU(2)) we see that the integrand depends on iV^^ only so that the 
integrals over dNf are trivial and give unity (upon proper normalisation). The partition function therefore reduceds 
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to 



z{K.) = / n '^atAa^t n d^^t/d^r/ n '^pf ^pff n n 

(tj) f t fct 



In 





X E n 

- E n 



t /ct 



1 

1 + - 

7 
1 

1 

7 



(3.34) 



where we also performed the integral over p ^ . 

Next we perform the integral over AN^^ to solve the simplicity constraint (implementing Eci.(3.22|) 



{a.t) 

X E n 

X E n 



z(/c) = /nd5^td5-,nd^*/ndp/(p/)'nd"*nd<^*/n'^(E/^/^*/^3^^^^^ 



1 

7 
1 



(3.35) 



where we also have performed the translation N^^ N^j:h^^f. Performing the translation — )■ g^-t'^t j the integrand 
no longer depends on and the Ut integral gives unity, leaving us with 



1 

1 + - 
7 



iT.t) (tj) f 

X E n 

X E n z^.,, 

Recall that for any SL(2,C) matrix g 



z{K) = / n '^atAa^t n dA^*/ n ^pf ^pff n d'/-*/ n e pi^tr^N^f' 

Pf) t^'j;^ {T'Nt'fh^^,^ec73e~^h^^^N;-/g~g-^,^ 



1 

1 

7 



(3.36) 



a b 



the representation matrix element 7r^„(g) reads 



V (j + (j - m)! (j + ?l)! (j - n)l ^j+n-i ^/ ^i-m-i 



(j + n- ly. {m-n + ly. {j - m - ly. l\ 



Applying this to ia^ 



10-3 = 



i 
~i 



(3.37) 



(3.38) 



(3.39) 



yields 



V U + my. {j - my. {j + n)! (j - n)! ^. 
(j + n - /)! (m - n + ?)! {j-m- ly l\ 



■j+n-l (^m~n+l qI ^_^y-m-l 



(3.40) 
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Likewise for 



■ 

e-'' 



we obtain 



^mn{f^4>) - Z^^j+n-iy{m-n + iyij-m-iyi\^^' " 
We conclude (summing over repeated indices), 



tr.+ [iNf, 



Since ie(j^e~^ = —ia^ we have similarly for the anti-self-dual part 

T^rnnl «^3j - (j- ^ _ ^ ^ _ ^ I 

thus 

We insert these formulae into the partition function 

^(^) = / n dff^-td^;-* n dp/ (p/)' n d^*/ n ^-^t/ n ^ PfN^f^N-^') 



1+ 



(tj) 



(tj) 



set 



Pf 



a,b,c 

and perform the integrals over d(f)tf which enforce = b'^j^ = b„f, and restrict the range of the sum over 
set • • • n {-j~f, ■■■ ,j~f}- Accordingly, 



Z{IC) 



n '^9:td9-t n dp/ (pff n ^^^z n -^(e /^/^'/^^iv-/ 

>,t) / ' 

En%,( 



(«,/) 

Pf) /5, 



t /ct 
1 



E\A ,(^''/)-f;,. (^r/)v^„. to, 

a,b^c 

':-\iN,f)^;' {N-^'W;J^_ {g^^g-,). 

"•af"<'f "<rl^<7f '-<T/"er/ 
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3.3 Topological/Gravitational Sector Duality, 7-Duality 

Before performing further computations, in this subsection we consider the topological sector I± of the simplicity 
constraint. Because we consider the model with finite Barbero-Immirzi parameter, the sector I± is actually also 
gravitational here in the following sense: By definition, tr(F A *(e A e)) is the Palatini (gravitational) term while 
tr(i^ A (e A e)) is the topological term. Since we are considering the Plebanski - Hoist Lagrangian tr(F /\{B + ^ * B)), 
inserting the gravitational solution _B = ± * (e A e) yields (due to *^ =id in Euclidian signature) the Palatini - Hoist 
Lagrangian with Immirzi parameter 7, that is, ±tr(i^ A (*(e A e) + A e)) while inserting the topological solution 
B = ±e A e yields Palatini - Hoist Lagrangian with Immirzi parameter I/7, that is ±itr(_F A (*(e A e) + 76 A e)). 



rescaled by I/7. If we change coordinates from to ±X^/"f in the partition function (2.1 1 then we obtain the 
relation 

Z^(/C) = /F-21T Zy^{IC) (3.48) 

where F, T respectively denote the number of triangles and tetrahedra respectively in JC (the powers arise from the 
Lebesgue measure and the 5 functions respectively). The appearing power of 7 drops out in correlators, hence up to 
the rescaling of the n-point functions of involving X^ , Zj, Zi/^ yield the same correlators. It follows that the model 



(2.1) is a mixture of gravitational and topological sectors as it should be. 

This is before restriction to either the gravitational or topological sector respectively and the manipulations (drop- 
ping infinite constants) that followed. For comparison, the partition function for the topological (I) and gravitational 
sector with Immirzi parameter 7 respectively read (before expanding the exponentials) 



Zi^^^'ilC) = f [U d^X+ d^Xj] [H dg+ dg-,] [J] dg+f dg^^] 
f i'y,t) (tj) 



/ct t (tJ) 

X exp(z[l + 7-I] J2 Tr(X+«;+^) + i[l - 7-I] ^ TiiXjuj-f)) (3.49) 

The only difference is the sign in the 5 distribution enforcing the linearised simplicity constraint. Now change variables 
X^ — > /7 in the model I (this induces also Xfj — >■ '^X^^ /j) . This switches the sign of the simplicity constraint 
to that of the model II, maps the I/7 in the exponent to 7 and rescales the Lebesgue measure and the S distributions 
according to 

Z^ifC)=-f^^-'^Zi}^{JC) (3.50) 

The power of 7 again drops out in correlators and thus up to 7 powers coming from n-point functions, "topological" 
correlators with respect to 7 are essentially the same as "gravitational" correlators with respect to I/7. We coin this 
relation between the two sectors "7 duality" . We will therefore not discuss model I any further in this article. 



4 New Spin-foam Model 

4.1 A Simplified Model without Closure Constraint 

In this subsection we discuss a simplified model by removing the closure constraint in the partition function Z{1C) 
by hand as it is also done in existing spin foam models. We do this just for a better comparison between our model 
and those models as far as the modifications are concerned that result from commuting rather than non commuting 
B fields. The discussion of the full model and the additional modifications that come from a proper treatment of the 
closure constraint will follow in the subsequent subsection. 
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The simplified partition function reads (from Eq.(3.47)) 



(a,t) 



1 - 



Pf 



a.h,c 



(4.1) 



In order to explore the structure of the spin-foam amplitude (e.g. vertex amplitude) for this partition function, we 
use the following recoupling relation {N £ SU(2)): 

= b+\(g>{r,b'\ U{N) \j+ , a+) ® ir , a-) 



J2 {k,a+ + a-\j + ,a+;r,a-) {j + ,b+ ; F ,b-\k, b+ + b-) tt^^ (N) 

k=\j+-j-\ 



(4.2) 



We denote by c(A:,j*)° " = c(fcij )"+a- Clebsch-Gordan coefficients {k,a\j^,a'^;j ,a ), which are real and 
vanish unless a — -\- . Thus (summing repeated indices) 

(iv)<:,, {N) 



3 ' +] 



J2 c{k,3^)a-^a^<k,J^)l,-<m 



(4.3) 



By using this recoupling relation wc find 



^ c(fc,J±)^+,_ TT^^C^t'/) c(fc>J^)M c(fc',J^)M 

fe,fc' = |j + -J-| 



(4.4) 



where (5 and a' are fixed to be 26. We note that k and fc' are restricted to be greater than or equal to \2h\ which 
we take care of by defining c(fc, j*)^^{, to be zero when k < \2h\. Inserting this result back into the partition function 
■Z'simpiificd(/C) resuhs in 

■^Simplified (A^) 







^F ( 






( 


1+1 


1-1 






7 




7 



i'^.t) f (tj) {i±.}('^,/) 

a} a-f (Tf cr/ erf erf erf a 



E _ 



(4.5) 



Now we consider a vertex v dual to a 4-simplex cr. We fix the orientation of each dual half edge (tr, t)(= (w,/x) in 
the notation of FIG. p| to be outgoing from the vertex and integrate the SU(2) holonomies gJt- The integration of 



leads to a result that depends on the orientations of the wedges bounded by (cr, t) . We say a wedge w bounded by 



{<j,t) is incoming to the edge (cr, t), if the orientation along its boundary agrees with (cr, i), otherwise we call it to be 
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Figure 2: A wedge (cr, /) with a S0(4) representation [j'^^^j^j) and two SU(2) representations fco-/ and fc^y, where 



outgoing from (ct, t). The integrations of g^^ in Eq.(4.5) 

y dg^t (g) (s^t) 



(4.6) 



w incoming {(J.t) 



w outgoing (cTji) 



equals a projection operator "^^t for each dual half edge (a, t 



OS) ^.i OS) 

It' incoming w outgoing {o',t) 



v*+ 



OS) OS) 

incoming (ct,^) outgoing (cr,t5 



Inv 



^.i OS) 

incoming {(T,tj w outgoing {(T,t\ 

t 



I Inv 



V- 



(8) 



It! incoming ((7,i) w outgoing {cr.t) j 

t 



(4.7) 



where Vj is the representation space for SU(2) unitary irreducible representation, and we keep in mind that each pair 

in 



(cr, /) determines a wedge and C^t (jct/) ^'^^ 4-valent SU(2) intertwiners forming an orthonormal basis ir 



Inv 



V*± 



(4.8) 



yii! incoming j w outgoing (cr.i) y 



Thus the resuh of the integrations of g^^ in Eq.(4.5 ) is a product of the projection operators ^p^t for aU the dual half 



edges (cr, i). According to the index structure appearing in Eq.(4.5 1, we find that in each *Po-t the adjoint intertwiners 



C'^± (jct/) ^'"^ combined with the indices a^^, c^^, where are for the incoming wedges and c^j are for the outgoing 
wedges. However the intertwiners C^'i (j^/) f°^' each half edge (cr, are contracted with other half edge intertwiners 
of (cr, t'j at the vertex dual to cr. Summing over the indices a^pC^j,, the integrations of g^^ in Eq.(4.5) result in a 
product of 



yc% c% ■ yet (j-,) (j-,) 



(4.9) 



for all half edges (cr, t), where • • • are the indices contracted with other half edge intertwiners of (cr, i') at the vertex 

;he intertw 



dual to cr. According to the structure of Eq.(4.5), we assign the intertwiners 
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Figure 3: v is the vertex dual to the 4-simplex a. /x = 0, 1, • • • ,4 are the labels for the tetrahedra forming the 
boundary of a. The edges [v, fi) fi = 0, ■ ■ ■ ,4 are the edges dual to the tetrahedron t^. The face determined by v, /i, i' 
fijV — 0, ■ ■ ■ , 4 is the wedge determined by a, t^^t^- 



to the beginning point of (cr, i), while we assign the adjoint intertwiners 



X^^f) {a+A.{c+^,} ' '^tt (■^'"■^){a- }.{c-,,} 



(4.11) 



to the end point of (cr, t). 



The contractions of the half edge intertwiners Eq.(4.10) at each vertex dual to cr gives a S0(4) 15j-symbol 



:= tr 





tr 








{cr,t) 



(4.12) 



to each 4-simplex a (to each vertex dual to a). 

On the other hand, each of the adjoint intertwiners Eq.(4.11 ) at the end point of (cr, t) is contracted with the Clebsch- 
Gordan coefficients c(/c, j'*)"+^_ and c{k' , . Thus we obtain a 4-valent SU(2) intertwiner {kcf,k'^j;j^^) at 



the end point of each half edge (cr, t) (summing repeated indices) 



<T,J) incoming {<t,j ) outgoing 



(4.13) 



If we choose an orthonormal basis in the space of 4-valent SU(2) intertwiners (labeled by lat) 

I 



OS) ^^^f Qs) 

incoming (cr.i^ /' outgoing {(Tjt] 



(4.14) 
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we may expand itt [kaf, k'^f',]^^) in terms of this basis, explicitly 



E 



, {P<T/},{p'„f/} 



Insert these findings into the partition function ^simplified (A^) yields 

■^Simplified (^) 

- E E/ndP/(P/f ndiv./n%, 



(4.15) 





1 


( 


1 




( 


1 + - 


1 - - 






7 




7 



?n{i«}3„«(. 



((T,t) 



n -iz.2,j^^'f) [<k.f,j%)Z%^ <k'^f,j^f)Z%, 



'26, 



from which we we read the vertex amplitude for each vertex dual to a 4-simplex a 



(4.16) 



(4.17) 



(CT,t) 



Next, we consider the integrations of dNtf. Since the closure constraint is removed, the integrals over dNtj can be 
done immediately. Consider a tetrahedron ti is shared by two 4-simplex a^, ai+i (c.f. FIG{T|, the integral of dNt-j is 
essentially 



There are three consequences from these integrals: 



dim(fc^,/) 



i+if 



(4.18) 



1. The SU(2) representations ka--f is identified with fc^.^^y, thus we label 

= fc^.+i/ = kt,f 

where ti is the tetrahedron shared by the 4-simplices ai,ai+i (see FIGB. 



(4.19) 



2. For the SU(2) intertwiners on the half edges {<Ti,ti) and ((Ti+i,t 



n 

/ incoming {<Ti ,(T^-|_i ^ 



n 



/' outgoing ((Ti,(Ti+i) 



which identify the half edge intertwiners into full edge intertwiners 



7 - =7 



(4.20) 



(4.21) 



where := (cTi, cr^+i) is the edge dual to the tetrahedron ti. 



26 



Figure 4: SU(2) representations ktf assigned to each pair of {t,f ). 



3. For each face dual to /, we have a factor 

kl/ 



kl/ 



kl/ 



bf i=i 

where the indices b^f are identified for the different wedges belonging to the same dual face and the ran] 
is 

kl/ 



1=1 



and \(7\f is the number of vertices around a face dual to /. 
Finally we consider the integrals of dp/. We define a triangle/face amplitude 



A 



■tf 



kl/ 



bf 4=1 



1 

kl/ 

dp/(p/f n 



i=l 







^0 ^^./ ( 


1 




.^^■:./ ( 


1+1 

7 


1 

7 





By Eq.(2.18), we can directly compute the expression of the function f3j 



y2,+i(2r) 



where Jn{x) is the Bessel fimction of the first kind. The proof of this relation uses the recurrence relation: 

J2j+2(2r) + J2j(2r) = (2^ + 1)' 



In the large-j limit j — >■ oo, the asymptotic behavior of the function Pj is 

Large-j : Pj{r) ~ r^^" ^^ 5{2r- 2j) 



(4.22) 
re of 5/ 

(4.23) 



(4.24) 



(4.25) 



(4.26) 



(4.27) 
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It follows that in the limit j^j -> oo, the asymptotic behaviom^ of the Bessel functions constrains the S0(4) represen- 
tations on the wedges by 



and also impose the well-known constraint on the self-dual and anti-self-dual representations 



1 



1 



(4.28) 



(4.29) 



which gives the quantization condition for the Barbero-Immirzi parameter 



If |7| > 1 : 
If |7| < 1 : 



7 



7 



^7 -J/ 



(4.30) 



While it is nice to see that we obtain certain points of contact with the EPRL and FK models respectively, one should 
keep in mind that these constraints hold only in the sense of large-j. In general, the representations which do not 



satisfy Eqs.(4.28) and (4.29) still have nontrivial contributions to the spin-foam amplitude and it is not clear whether 
these "non EPRL/FK configurations" have large or low measure. 

Let us summarize the structure of the partition function ^simplified (^) in terms of the 4-simplex/vertex amplitude, 
tetrahedron/edge amplitude and triangle/face amplitude 



^Simplified (/c) = ^ E E n ^/ (^'^z' ^*/) n Mhf) n {jf/^^fje 

{j±jR,}{ie} / t 

where kt / is constrained by the condition that for a tetrahedron t shared by both 4-simplices a, a' we have 



ktf e 



(4.31) 



(4.32) 



and the 4-simplex/vertex amplitudes, tetrahedron/edge amplitudes and triangle/face amplitudes are respectively iden- 
tified as 



Mhf) 



n ^ 



(a,t) 



Mi 



bf i=l 



Pf ^ 



Pf 



When we take the large-j limit: j^r — > oo, by the previous discussion, we obtain the constraints; 



Jaf = Ja'f = jf and 



1 - 



J/ 



(4.33) 



(4.34) 



Thus in the large-j limit the spins j^y for different wedges are identical on the same face dual to /, and and 
satisfies the "7-simple" relation in this limit. Then the vertex amplitude reduces to 



{it} 



(4.35) 
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where and jj subject the relation in Eq.(4.40). We notice that in the large-j hmit Eq.(4.35) is nothing but 
the vertex amphtude of the FK^ Model (when I7I > 1) [TT]. And in the large-j limit the integral over pf can be 
approximated by a discrete sum over jJ or in the path integral Eq.(2.1 1. 



4.2 On the Implementation of Closure Constraint 

In this subsection we properly keep the closure constraint in the partition function: 



{3 



n <Z^2KAN^'f) [<k.f,j^f)Z%, c{Kf,j^f)Z;Q 4:„,;,(^r/) 



Pf 



(4.36) 



Here we can also extract the vertex/4-simplex amplitude Aa-, the edge/tetrahedron amplitude At, and the face/triangle 
amplitude Af 



/ nd^*/'^(E'°/^*/^3A^r/) 



/ct 



/ct 



/ outgoing /' incoming 



1 

1+ - 
7 



1 

1 — 
7 



Pf 



Then the partition function can be written in terms of these amplitudes as: 

•^f [pf'^iaf-'^'yf-'Kf'-.baf^ At(^pf]kaf,Kf;huh„f^ Aa (^j^j:;k^f,k'^f;lo 
Note that in the large-j limit, 

^ .4- 

1 



(4.37) 



Large-j : /3.±^ 









( 








7 



1± 



P/ 



(5 2 



1 ± 



7 



(4.38) 



(4.39) 



In this limit, the integral of pf is completely constrained by the delta functions. Thus, as in the previous section, the 
delta functions impose the constraints: 



jaf = Ja' f = jf and 



1 



1 



^7 



1 



1 



J/ 



(4.40) 



29 



This shows that in the large-j Umit the spins j^j, for different wedges are identical on the same face dual to /, and j'j 
and jj satisfy the "7-simple" relation in this limit. However, since at the current stage the constraint 

k.a = K.+,f ^ Kf (4.41) 

are not obviously imposed by the integral of Ntf (because of the present of closure constraint in At), the vertex 
amplitude Aa, even in the largc-j limit, does not approximate the FK-y vertex amplitude in general. 

To explore the structure of this amplitude, we consider the integral of Nf f in the expression of At, for a tetrahedron 
t shared by a, a' 

■^t (yPf ; kaf , Kf ; i^t ; &<7/) 

:= I W^Ntf 5(y" pfNtfTiN-A Cf (k^f,k',;) Cf (k„'f,k'f,) 

n 4;„../A^r/)-fc;,2^,,w/) n -ar;,2^,,w/')4:f;,,.,,,(^t7^) (4.42) 

/ outgoing /' incoming 

Recall that Ntf = gtfNf {Ntf = N^j = g^j:N^) while the integrand of the partition function only depends on the 
combination NfT^Ny^ (recall that the integrand depends on X^), i.e. the integrand is invariant under Nf i— )■ Nfh^ 
where € U(l) thus it only depends on SU(2)/U(1). Let us parameterize Nf in terms of the spherical coordinates: 
In terms of the complex coordinates {zf, Zf) on the unit sphere we have 

N{zf) = -jf^i \ "l] (4.43) 

where the complex coordinates z, z are defined by the stereographic projection, and the unit vector on S"^ is expressed 
in terms of the complex coordinates 

*) = -■(-^n^^<''+fi^''^+^<"')='^w^«'^w- ("^) 

Under the action of SU(2) group 



where 



Therefore 



9N{z) = N{z^) ^ where I--' ^ e U(l) (4.45) 

\ " \a-bz\ J \ ^ \a-bz\ J 



az + b . , I a b \ , . 

with g=\ _ _ (4.46) 



a-bz \ -I) a 



Ntf = 9tfN{zf) = N{ztf)h-^l^ V, e U(l) (4.47) 

where ztf = ^J*^. Note that the above decomposition may also be understood by writing the SU(2) matrix in terms 
of Euler coordinates, i.e. u = u{4>2) u{9) u((/'i) for all u & SU(2) 

where u{(j)\), u{<p2) € U(l), and < ^i, (?i2 < 27r, < ^ < tt, while the SU(2) Haar measure can also be written as 

dg = — sin0d0id6'd02 (4.49) 
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Hence the integral Eq. ( 4.42 ) can be written as 



Af(^Pf'iKf,Kf'-:^<Tt;Kf^ 

/ outgoing 



},{/3„'/}, 



(4.50) 



/' incoming 



where AD.tf is the standard spherical measure on S"^ = SU(2)/U(1). Since i^lan{^<t>) = (e 



2im4 



it follows 



/ct 



/ outgoing 
/' incoming 



(4.51) 



The integrals /(^^^^ impose the constraint that b^f — b^if = bj for all f dt, hence 



kcr'f,k'^,f, 



/ct 



/ct 



{^.'/'},{/9„'/} 



n 4;,.„,(^(-t/)-^) -fc;^2.,(^^(-t/)) n <;.2v(^(^*/')) 4;v,,(^(-*/')-0(4.52) 

/ outgoing f incoming 

Moreover for the outgoing dual face /, we have the relation 

- (-1)-^/-— -■/x^,^_,,^(iv) 4z,,,{n) 
= (-1)-^/-—-^ "e" ^(^^/;fc:^/'fc-/)-l,,„,/('^/;fc^/'^-/)-2.,2.,-;^;oi 

while for the incoming dual face /' we have similarly 



N) 



(4.53) 



= (-1)^'=-/'-^"/'-^-/' E ^ ^(W';^'^/''^:.'/'):':;,-,,,, K'*/';fc./'-fc;'/')2.,„-2v'^;:'^oW4.54) 

Thus the integral reduces to 

A{pf]k„f,k'^f,l^t]b„f^ = J Y[d'^^tfd(j)tf Sl^'^pf^tf^ et{latJa't,k'^f,K,f,Kf',k'^,f,,bf;ztf) (4.55) 



/ct 



/ct 
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with the integrand {t is the tetrahedron shared by a, a') 



f outgoing 



f incoming 



't/'=|fe<T/'-fe^/j,/l 



(4.56) 



The complete integration of Eq.(4.55l turns out to be rather involved, thus is left as a future research. What one 



can say is the following: From the expression Eq.(4.56), it is not hard to see that the set of amplitudes contributing 
to the simplified model .^simplified (^) are part of those of the full model Z(/C). To sec this, notice that the simplest 
nontrivial contribution of the integral in At comes from the term with Itf = — 0. With the constraints = l^j = Q 
(and thus p^y = ptf — 0), we obtain the same set of constraints as it was in the previous subsection for the simplified 
model ^Simplified (^)- 



Since 



( — 1)''^" 

-\/dim(fc) 

we obtain, by extracting the term with It /' = l^^ — and dropping the contribution from the other terms 

6t {lcrt,la't,Kfyf^c'f^f^crf':K'f'^bf; Ztf) 



(4.57) 



(4.58) 



f outgoing 



(-l)'=^/+«-/(-l)'=^/+2''/ 



/' incoming 



fc„,,-26. 



dim(A;o./') 



n dim(fct/)'^''-^''''-'^'^''''^' 



J at 

For this subset of amplitude the edge/tetrahedron amplitude reduces to 



(4.59) 



ndim(fct/)'^'^^/ 



/ct 



/ct 



(4.60) 



Then we can define a spin-foam model by pick out a subset of amplitudes in the full partition function Z{IC): 



(4.61) 
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The amplitudes in Z'{K,) arc contributions with the closure constraint implemented, however unfortunately they may 
not exhaust all the contributions. 



tween the 



In Eq.(4.60) the Kroneckcr deltas 5k' k ,,5^ k' 5, a imply that there is an one-to-one correspondence be- 
transition channels in the simplified model 2^simpiificd(^) f^nd the transition channels in the model Z'{IC), 
which form a subset of the transition channels in Z(/C). Consider the sets {^simplified} and {Z} respectively, which 
are the collections of spin- foams that contribute to their respective partition functions ^simplified (A^) and Z(JC). Our 
above analysis then reveals 

{^Simplified} C {Z} (4.62) 

At this point this is all we can say about the relation between the models with the closure constraint in place or 
not. The additional weights and contributions in the full model may severely change the correlators (physical inner 
product) and it is by no means obvious that the simplified model is a good approximation. 

As a final remark, the above inclusion is in terms of spin-foam amplitude, in the sense that we write the partition 
functions as a sum of amplitude over possible spins and intertwiners. Moreover such an inclusion is natural from the 
path integral point of view. We consider a simple example: Consider a function f{x, y) on which has a Fourier 
transform f(k, q) and that we have a "closure constraint" y — Q. Then the Z integral (with closure) corresponds to 
(dropping factors of 27r) 



Z = j dxdy5{y,Q)f{x,y) = j dxdyS{y,0) J dkdqf{k,q) ejip{i{kx + qy)) 
dx J dkdqf{k,q) exp{ikx) = J dkdqf{k, q)S{k,0) 

dqf{0,q) (4.63) 
On the other hand the ^simpiied integral without closure is 

•^Simpiied = J dxdyf{x,y) = J dxdy J dkdqf{k,q) exp{i{kx + qyj) = J dkdqf{k,q)6{k,0)S{q,0) 

= 7(0,0) (4.64) 

Hence the Z amplitudes are more in Fourier space {k,p) corresponding to spin-foam representation, and less in real 
space (x, y). 



5 Outlook 

In section|4]we first carried out the analysis for the simplified partition function without closure constraint and obtained 
the spin-foam model .^simplified (^)i then we discussed the complete partition function Z(^K.) with closure constraint 
implemented, however we did not compute yet explicitly the full set of possible spin-foam amplitudes. We were only 
able to show that all the spin- foam amplituded contributing to ^simplified (A^) are contained in those contributing to the 
full model Z{)C). Therefore, in addition to present spin foam models, our commutative B field model variable sums 
over additional amplitudes having non- trivial contributions to the partition function Z{IC). While we have shown 
that in the large-j limit the 4-simplex/vertex amplitude of .^simplified (^) can be related to the 4-simplex/ vertex 
amplitude of FK-^ Model (I7I > 1), for the full model Z{IC), even in the large-j limit, there exist additional, non-trivial 
spin-foam amplitudes. It would be important to further specify those unknown spin-foams contributing to {Z} but 
not to {^Simplified}, at Icast for their large-j asymptotics. 

Unfortunately, the relation between our new model and EPRL model is almost untouched in the present article. 
Although we have seen that all the EPRL spin- foams (with possibly different triangle/face and tetrahedron/edge 



33 



amplitudes) are included in {^^simpiificd} (thus in {^}), it seems to us that, however, they are not quite special among 
the spin-foam amplitude contributing ^simplified (A^) or Z{IC). We expected that the relation between our model Z{IC) 
and EPRL Model could be realized by the non-commutative deformation, like in the case of Barrett-Crane Model. 
The reason for our expectation was that (1) both models are defined via the non- commutative operator constraint 
technique, and (2) when the Barbero-Immirzi parameter 7 — > 00, EPRL Model reduces to Barrett-Crane Model. 
However it turns out that our expectation is difficult to realize, since the non-commutative deformation via the group 
Fourier transformation hardly works for the case of finite 7. It seems to us that if our model Z{IC) and the EPRL 
Model could be related via any non-commutative deformation, we should rather choose a different deformation scheme. 

The present article starts from a purely path- integral/spin- foam point of view. If we also consider the relation 
between the path integral and canonical quantization, then the partition function Eq. (2.1 1 should probably be 
modified. It is pointed in 13^ that a quantum gravity path integral formula consistent with canonical physical inner 
product should not only be an naive path integral Eq.(1.2) of Plebanski-Holst action, but also include a suitable local 
measure factor in the path integral formula. The local measure factor is a product of a certain power of spacetime 
volume elements and a certain power of spatial volume elements at all the spacetime points. The implementation 
of such local measure factor in the partition function will modify both the 4-simplex/ vertex and tetrahedron/edge 
amplitudes. A detailed analysis of this issue will be postponed to future research. 
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A Noncommutative Deformation and Barrett-Crane Model 

A.l Noncommutative Deformation 

In order to further investigate the question, in which sense the closure constraint is redundant when working with 
non commutative B fields (as is common practice in existent spin foam models), in this section, we explore a non- 
commutative deformation of our starting point, the partition function Z{]C) in Eq.( |2.1| . The non-commutative defor- 
mation we will employ here comes from a generalized Fourier transformation defined on a compact group |24) . The 
deformation replaces the normal c-number product in the expression of Z{IC) by a non-commutative "★-product" (we 
will briefly review the definition below). Surprisingly, this non-commutative deformation establishes a relation between 
the new spin- foam model Z{IC) we analyzed in the previous section and the Barrett-Crane spin-foam model [9]. 

First of all, we recall the partition function Z{IC) in the commutative context (after the linearization of simplicity 
constraint): 

Z{IC) := / n d'^;d^^7 n d5+td5-, [J ^5^/^5*7 H H ^ (^t/ + P^tX+u^') I] ^ ( E ^*/) 
/ {<y:t) (tJ) t tj t fct 

X Yl e'*'(-^/^'/'S''-»^t'»?/) Jl g«f(^797t9t.ff;t'9;/) (A.l) 

where (3 — and for the convenience of the following analysis, we have made a change of variables 

Xf^{l±]^ Xf. (A.2) 
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and dropped a constant 7 dependent factor. Here we assume that our structure group is SO(3)xSO(3) instead of 
SO(4)~SU(2)xSU(2)/Z2. The reason for this replacement is to be compatible with the group fourier transformation, 
which will be seen shortly. 



We now replace (by hand) the commutative c-number product in Eq.(2.1) by the non- commutative ★-product on 
su(2) ~ M3 defined in 'M], that is 

^4^t^(X\gi\) ^ g5^tr(X|g2|) pJstr{X\gig2\) ^^ 3^ 

where a is the deformation parameter, X — X^Tj and tj = —icj with aj the Pauli matrices (Ticrj = Sij + ieijk<^k, 
g G SU(2) represented by a 2 x 2 matrix and \g\ — sgn{ti-g)g so that \ — g\ ~ \g\. We can write g G SU(2) as 

g = P(i + iaP -a, Pq^ ^2| |p| |2 ^ ^-^^^ 

Thus l^l is the projection of g on the upper "hemisphere" of SU(2) with Pq > 0. Therefore the "plane wave" in 
Eq.(A.3) can be written 

depends on SO (3) only (its character expansio depnds on integral representations only because it is an even function 
under reflection g -H- —g) . With these "plane waves" we can define an invertible "Group Fourier Transformation" from 
the functions f{g) on S0(3) {f{g) — f{—g) for g G SU (2)) to the functions f{X) on the Lie algebra su(2) 

f{X) = f dg fig) egiX) 



fig) = j d?X f{X)^eg-.{X) = ^ — j d'X fix) eg-.{X) (A.6) 

Given two functions fiiX) and f2iX) in the image of the group Fourier transformation, their ^-product is defined as 

/i(X)*/2(X) = j d5id.g2 /i(gi) /2(52) e,,, (X) * e,,, (X) (A.7) 

and when the deformation parameter turns to a — > 0, the ^-product reproduces the normal commutative product (if 
we keep Pq.P fixed, see (A.4)). 

We also have two identities for delta functions 

^so(3)(<?) - 8^ / <^'Xe,iX) 

5xiX') = J dg eg-^iX) egiX') (A.8) 

where the second delta function is the Dirac distribution in the noncommutative sense, that is 

J d'X' iSx * /) iX') = J d'X' if * Sx) iX') = fix) (A.9) 

With the above definitions, we can make a noncommutative deformation of the integrand in Eq.(A.l). In the 
following we fix the deformation parameter to 

a = £l = l (A.IO) 
The reason for this choice is that only in this case the closure constraint turns out to be redundant and can be removed 



from Eq.(A.l), which is necessary in order to derive the Barrett-Crane model. We will show this immediately in the 
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next paragraph. On the other hand, fixing a = makes it impossible to study the commutative limit a — )• of the 
non commutative model Za{IC) which we denote by and thus we cannot compare with the commutative model Z{IC). 
We first define the noncommutative deformation of 

n<E^*/) n e^*<"^'^-^-^-'^''^) = n<E5*!/^/V*) n e^"("^'^-^-^-^^^) (a.ii) 



* fCt (ct,/) t fct (aj) 

Given a face dual to the triangle / with n vertices dual to the 4-simplices cti, • • • , ct„ (cf. FIGjl]), we define the quantity 

(A.12) 



*---*e+ ^ + -k e + + + + 



A possible noncommutative deformation of Eq.(A.ll) is 



l[^htl[g+iX+,g+,g+,ht) 



t f 

because the noncommutative Dirac distribution for the closure constraint is 



/ct /ct 



(A.13) 



(A.14) 



It is here where the choice a ^ Pp was important because we have implicitly set = 1 in the exponential so far 
(it comes from the fact that the flux field has dimension cm^ and the Plebanski action is multiplied by 1/k where 
kH = i^p) so restoring it we can combine the ordinary product of exponentials into star products only if the deformation 
parameter is given by a = £p. 
However, since 



g+{X+,gt„gtf.ht 



h n+ n+ n+ * ^ n+ h n+ n+ n+ * ' ' ' * 6^+ ^ ^+ ^+ „ + 



we can absorb ht into gta by a change of variables 



{X+) (A. 15) 



(A.16) 



while dgjr^ does not change since it is Haar measure. Therefore finally the integral of ht gives unity, which shows the 
redundancy of the closure constraint for this particular non commutative deformation! 
Next we consider the simplicity constraint 



5{X;j + putX+Ut^) = 5{g-fXJg^^ + /3ut5+ ^+5+ u,-') 
whose noncommutative version is 

5{gtfXjg-^ + ^Utg+X+g+Ut^) = 



(A.17) 



dvtf e^^jig^jX^ g^^ + (3utgtfX+ g+^u^ ^) 
•ivtf ey^f{g;'j:XygJ^) ey^j{Putg+fX+ g+^u-^) 



(A.18) 



where 



e^AX) := egipX) 



(A.19) 
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Thus for each face dual to the triangle / with n vertices dual to the 4-simplices cti, • • • , tT„, we define 



I'^{X+,g+t,g+,ht,Ut,Vtf,(3) 



+ , + *e++ + + * e J. 



t2 "'2 



*e+. + *e+ + + + *e I _i + * 



and 



-k e + , + -k e + + + + *e + _i . 



.^9/t,.9t„<,i9„ltl94j/ * %ti^'tl/ffti/ *%ti9ti<.2f-2t2ft2/ *%*2'"'2/9t2/ ^ 



(x;()\.2o) 



(x-f 



(A.21) 



Then the deformed partition function is defined by 

^4^) =111 d'^/ n d.9i^(d5-, n d5tV% n d/^^dwi n d^^*/ 



f 



(A.22) 



which is the noncommutative deformation of Eq.(A.l). However since we have shown the redundancy of the closure 
constraint in Zi,{IC), we can equivalently write 



uic) ■■= / n d'^;d3x7 n d.9+ dg-, n '^at^j n n ^^^z 

n -^f ^9tt^9tf.uuvtf,f3) FJ (Xy , 5" ,gyj:,vtf) 
f 



(A.23) 



where J^j" is replaced by 



J'f{Xf,g^t,gtf,ut,vtf,l3) 



e„+ „+ „+ „+ * e 



9ft„9t„aia„^tigtif aft^-^ti ^t^futiat^f 

■*:€++ + + *e'^i_i I 

9ft^ati^29<'2t29t2i 9/t2"t2 '"*2/"*29t2/ 

★ • • • ★ 



* e„+ 



9ft„_l9t^_-^a„aantn9i„f " " a^^^t^^inf^tnat^f 



ixf). 



(A.24) 



A. 2 7 = oo and Barrett-Crane Model 



The computation with general /3 is difficult, because it involves the T^-product between two different types of plane 
waves eg and e^, which is even not well-defined in general (since they could be consider having different deformation 
parameter). Therefore here we only consider the simplified case that 7 = oo. Then 13—1 and in this case we can 
directly compute Tj' to be 

^^{X^,9tu9tf'^t-,vtf,l3^ 1) 

~ ^9"jt„9t,,^l9t^t^'^tl^tlJUtiat^a29t2t2'^72^t2S^t2---9t^_^a„9t„t,,^7^-"t„fUt,,gt^f^-^f ) 



S^„„iS^iti«t/^'tl/"n9t%2^'^2*2"t2^'"*2/"t2---9i_l<T„S^„t„"t„^"t>i/'"t. 



(A.25) 



Similarly for the anti-self-dual part 

^J{XJ ,9^t^gtpVtf) 



"9/t„9t„^i9<Tifi"ti/£'ti,j2f<'2t2'"'2/---9t„_i^„9<T„t„''tn/9t„/ 



9t^cr^9CT^t2^^tl/9t^cr2 9(72*2 



(^7 



(A.26) 
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We define the following changes of the variables 

9t,u ^ dt.uUu 9ua,+, ^ '^u^9ua,+^ ^ 9ft^^tnX+Ut^g+^ XJ ^ 9jtn^j9t„f (A.27) 

where for each face dual to /, a unique tn{f) is chosen as the base point of the dual face. Thus the partition function 
can be written as 



UlC) := Jl[d'X+d'Xjl[dg+,dg-,lldvt 



f {tj) 

/ 
/ 

We perform the integrals over X^, XJ and obtain 

Z4I^) = [U d9ttd9-t n 

n ^ {9tTi9t^t^VHf9ta29t2t2'"t2f ■ ■ ■ 9t-ian9t„t„Vtnf) 

f 

n ^ (5r„<Ti5aifiWti/ffria2ff<^2t2^t2/ • • • 9t„_,ar,ga„t„Ut„f) (A.29) 

/ 

Decomposing the delta function into the charaters (here j's are integers since g G S0(3)) 



%) = ^dim(j)tr,(ff) (A.30) 



yields 



z.(/C) := E / n '^a^tdg^t n d^*/ 

j± (^,0 («,/) 

ndim(j+) 7rli^+ U+.,5^,tJ i^tif)---\i (fltt-K^.flXtn) <+6+ (^*n/) 
/ 

n d™(j7) (9ir.a^9a,t^) • • • (3t"„_ia„9;„t„) ^'-f,- (^t„/) (A.31) 



/ 

We immediately extract the integral 



= j dl^t/ T'lia^ (V/) ^I-b- (-f^t/) = / d-t^f/ '^6ia+ («^?f«~') ^'j-b- M 

dim(jj) dimOj) 



where ^ o)'^^~^^~ ^^^^ ^ — 9 ^ ^"^^ ^''i'^) is a real matrix. We insert this result into the 
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partition function (ny is the number of vertices around a face dual to /) 
Finally we obtain the result 

Jf (<T,t) / i=l ^ 

which gives Barrett-Crane vertex amplitude [9]. This result is consistent with the work done by colleagues using 
completely independent methods [25] . 
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